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Abstract 

Before  vehicle  mobility  in  snow  can  be  reliably  predicted,  a  complete  under¬ 
standing  of  motion  resistance  in  snow  is  required .  This  report  examines  several 
aspects  of  wheeled  vehicle  motion  resistance  using  results  obtained  with  the 
CRREL  instrumented  vehicle.  Resistances  of  leading  and  trailing  tires  are 
examined.  Limited  data  are  presented  for  undercarriage  drag,  and  third  and 
fourth  wheel  passes  in  the  same  rut  are  initially  analyzed,  as  is  how  snow 
deforms  around  a  wheel.  For  the  CRREL  instrumented  vehicle,  a  trailing  tire  has 
a  resistance  coefficient  of  about  0.0 17  tor  snow  depths  less  than  about  22  cm. 
For  deeper  snow,  the  disruption  of  the  snowpack  caused  by  a  preceding  wheel 
causes  snow  to  tall  into  the  rut,  resulting  in  higher  trailing  tire  coefficients.  For 
larger  vehicles,  which  in  some  cases  have  trailing  tires  carrying  larger  loads 
than  preceding  tires,  the  trailing  tire  coefficients  are  on  the  order  of  0.048  and 
0.025  for  second  and  third  trailing  wheels  respectively.  Since  there  are  no 
trailing  tire  data  available  tor  these  larger  vehicles,  these  values  are  based  on 
nonlinear  regression  analysis,  which  includes  a  prediction  of  the  leading  tire 
resistance.  The  results  and  observations  of  this  study  are  applied  in  a  reanalysis 
of  the  towed  resistance  data  obtained  during  the  U.S.  Army's  Wheels  vs.  Tracks 
study.  An  improved  algorithm  is  presented  for  predicting  wheeled  vehicle  motion 
resistance  caused  by  snow. 


Coven  CRREL  instrumented  research 
vehicle  in  deep  snow. 
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PREFACE 

This  report  was  prepared  by  Paul  W.  Richmorid,  Mechanical  Engineer,  of  the 
Applied  Research  Division,  Research  and  Engineering  Directorate,  U.S.  Army  Cold 
Regions  Research  and  Engineering  Laboratory.  Funding  for  this  research  was 
provided  by  DA  Project  4A762784AT42,  Design,  Construction  and  Operations  Tech¬ 
nology  for  Cold  Regions;  Task  CS;  Work  Unit  M04,  Vehicle  Snow  Mechanics  (Shallow 
Snow). 

Technical  review  of  this  report  was  provided  by  Sally  A.  Shoop  and  George  L. 
Blaisdell,  both  of  CRREL.  The  author  expresses  his  appreciation  to  Byron  Yoimg, 
who  assisted  in  conducting  the  tests  and  in  the  operation  and  maintenance  of  the 
CRREL  Instrumented  Research  Vehicle. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional 
purposes.  The  citation  of  brand  names  does  not  constitute  an  official  endorsement 
or  approval  of  the  use  of  such  commercial  products. 
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Motion  Resistance  of  Wheeled  Vehicles  in  Snow 


PAUL  W.  RICHMOND 


INTRODUCTION 

Models  of  terrain-vehicle  interaction  require 
algorithms  that  can  predict  or  estimate  the  resis¬ 
tance  to  motion  caused  by  the  terrain.  As  a  vehicle 
moves  through  snow,  it  must  overcome  the  resis¬ 
tance  caused  by  snow  deformation  to  maintain 
forward  motion.  There  is  resistance  at  each  wheel 
or  track.  As  a  vehicle  moves  through  undisturbed 
snow,  the  front  wheels  move  and  compact  this 
snow.  Succeeding  or  trailing  wheels  move  over  or 
through  the  snow  compacted  by  the  front  wheels, 
presumably  encountering  a  smaller  resistance 
(since  the  snow  has  been  compacted  and  moved 
out  of  the  way).  If  the  vehicle's  undercarriage 
drags,  this  adds  another  source  of  resistance  that 
must  also  be  considered. 

The  motion  resistance  of  a  vehicle  moving 
through  snow  is  a  very  complex  phenomenon, 
studied  since  the  early  1940's.  However,  a  robust, 
simple  model  of  motion  resistance  due  to  snow 
has  eluded  researchers. 

Currently,  the  U.S.  Army  uses  in  its  mobility 
model  an  empirical  equation  for  motion  resistance 
in  snow  developed  by  Richmond  et  al.  (1990).  This 
equation  is  applicable  to  a  large  range  of  vehicles, 
both  wheeled  and  tracked,  and  does  not  require 
elaborate  measurements  of  snow  properties.  How¬ 
ever,  it  is  currently  the  weakest  link  of  the  overall 
Cold  Regions  Mobility  Model  (Richmond  et  al.,  in 
prep).  When  developing  this  equation,  we  ana¬ 
lyzed  a  large  quantity  of  total  vehicle  resistance 
data  by  dividing  the  whole  vehicle  resistance  by 
the  number  of  wheels,  thus  making  the  assump¬ 
tion  that  the  trailing  tire  resistance  was  equal  to  the 
resistance  of  preceding  wheels.  This  assumption 
was  based  on  attempts  to  correlate  data  in  which 
motion  resistance  was  measured  on  one  axle 
(CRREL  Instrumented  Vehicle  [CIV]),  on  two  ax¬ 
les  (High  Mobility  Multipurpose  Wheeled  Vehicle 
[HMMWV]),  and  on  four  axles  (Heavy  Expanded 
Mobility  Tactical  Truck  and  Light  Armored  Ve¬ 
hicle  [HEMTT  and  LAV])  during  towing  tests. 
This  approach  was  not  without  basis,  as  the  den¬ 


sity  of  the  deformed  snow  behind  each  succeeding 
axle  generally  increased,  although  axle  load  did 
not.  Density  measurements  in  the  wheel  ruts  (Table 
1)  were  always  well  below  what  is  considered  the 
critical  density  of  snow  (0.5  to  0.55  Mg/m^)  (Yong 
and  Fugue  1977),  the  point  at  which  no  further 
compaction  by  typically  loaded  vehicle  traffic  is 
observed.  Additionally,  Blaisdell  (1987)  made  some 
limited  indirect  measurements  of  trailing  tire  re¬ 
sistance  using  the  CIV  and  found  that  trailing  tire 
resistance  in  snow  could  be  as  high  as  84%  of  the 
leading  tire  resistance. 

In  an  effort  to  increase  the  accuracy  of  the 
motion  resistance  model  for  undisturbed  snow,  a 
field  study  was  conducted  using  the  CIV  with  load 
cells  installed  on  all  four  wheels.  Prior  to  this  work, 
only  one  set  of  load  cells  was  available  and  these 
were  installed  on  the  front  axle  of  the  CIV  (from 
1981  to  1991).  By  use  of  the  data  acquired  from 
these  and  other  tests,  a  further  understanding  of 
the  motion  resistance  of  wheeled  vehicles  in  snow 
was  obtained,  allowing  an  improved  prediction 
algorithm  to  be  presented. 

EXPERIMENTAL  PROCEDURE 

The  experimental  procedure  was  designed 
around  the  capabilities  of  the  CIV  (Fig.  1  [Blaisdell 
1983,  Shoop  1992]).  The  instrumentation  of  the 
CrV  was  recently  enhanced  through  the  addition 
of  load  cells  on  the  rear  axle  and  a  new  data 
acquisition  system  based  on  a  personal  computer 
(Shoop  et  al.  1991).  These  triaxial  load  cells  mea¬ 
sure  longitudinal,  vertical  and  lateral  forces  at  the 
tire/terrain  interface.  Longitudinal  forces  repre¬ 
sent  net  traction  or  motion  resistance,  depending 
on  whether  the  tire  is  driven  or  free  rolling.  Indi¬ 
vidual  wheel  and  vehicle  speeds  are  measured 
simultaneously  with  the  wheel  forces.  By  towing 
the  CrV,  resistance  forces  could  be  measured  at  all 
four  wheels,  and  a  tow  hitch  that  could  be  placed 
off-center  (Fig.  2)  allowed  the  CIV  to  be  towed 
through  undisturbed  snow. 
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Table  1.  Vehicle  sinkage  and  snow  density  data  (after  Richmond  et 
al.  1990). 


Vehicle 

Contact 

pressure* 

(kPa) 

Snow 

depth 

(cm) 

Sinkage 

(cm) 

Initial 

density 

(Mglm^) 

Rut 

density 

(Mgfm^) 

HMMWV+ 

101.8 

15.0 

8.0 

0.20 

0.38 

123.3 

15.0 

8.0 

0.20 

0.48 

101.8 

19.5 

11.9 

0.20 

0.47 

123.3 

19.5 

13.8 

0.20 

0.50 

101.8 

23.5 

19.1 

0.20 

0.535 

123.3 

23.5 

21.0 

0.20 

101.8 

19.5 

14.5 

0.20 

0.455 

123.3 

17.5 

13.1 

0.20 

0.470 

101.8 

18.0 

14.5 

0.20 

0.445 

123.3 

18.0 

13.0 

0.20 

0.475 

101.8 

13.0 

9.0 

0.12 

0.320 

123.3 

12.0 

8.5 

0.12 

0.410 

101.8 

19.5 

14.0 

0.12 

— 

123.3 

18.0 

14.0 

0.12 

— 

123.3 

10.0 

9.0 

0.19 

0.575 

HEMTT* 

233 

20.0 

12.4 

0.25 

0.490 

230 

21.0 

16.9 

0.25 

0.440 

225 

20.0 

16.2 

0.25 

0.520 

218 

20.0 

15.6 

0.25 

0.460 

233 

16.5 

12.4 

0.25 

0.490 

230 

18.0 

13.2 

0.25 

0.490 

225 

20.5 

15.4 

0.25 

0.510 

218 

20.0 

15.6 

0.25 

0.510 

233 

20.0 

14.3 

0.245 

— 

218 

17.0 

12.6 

0.245 

0.460 

233 

13.0 

10.0 

0.12 

— 

HEMTT+t 

177 

17.0 

12.6 

0.245 

0.510 

174 

18.0 

15.1 

0.245 

0.525 

224 

17.5 

14.6 

0.245 

0.430 

212 

16.0 

13.1 

0.245 

0.455 

212 

16.0 

12.2 

0.245 

0.455 

212 

4.0 

3.0 

0.095 

— 

177 

12.0 

10.5 

0.12 

0.300 

212 

11.0 

9.0 

0.12 

0.315 

*Contact  pressure  of  the  wheel  behind  which  the  rut  density  measurement  was  made. 
tWith  an  inflation  pressure  of  137.9  kPa  (front)  and  151.7  kPa  (rear). 

"With  an  inflation  pressure  of  241.3  kPa  (front  two  axles),  275.8  kPa  (rear  two  axles). 
t+With  an  inflation  pressure  of  137.9  kPa  (front  two  axles)  and  206.8  kPa  (rear  two  axles). 


Tests  were  conducted  at  the  Keweenaw  Re¬ 
search  Center  of  the  Michigan  Technological  Uni¬ 
versity,  in  northern  Michigan,  during  January  and 
February  1993.  The  testing  areas  had  been  graded 
free  of  vegetation  prior  to  any  snowfall,  and  were 
nearly  level.  A  layer  of  packed  snow  was  beneath 
the  undisturbed  snow  in  which  the  tests  were 
conducted.  The  CIV  was  equipped  with  Michelin 
LT  235XCH4  tires  for  diese  tests. 

Initial  tests  were  done  with  the  tires  inflated  at 
103  kPa  (15  lb/in.2),  while  most  of  the  tests  were 
run  at  an  inflation  pressure  of  179  kPa  (26  Ib/in.^). 
Tire  data  are  in  Table  2. 


A  static  calibration  procedure  (Shoop  1992)  was 
used.  Prior  to  coUecting  data  in  the  snow,  we 
conducted  at  least  four  tests  on  hard  surfaces. 
These  hard  surfaces  were  sometimes  covered  with 
packed  snow  or  ice;  however,  an  attempt  was 
made  to  use  areas  that  were  as  smooth  as  possible. 
These  hard  surface  data  were  averaged  for  each 
load  cell  and  used  to  remove  from  the  subsequent 
snow  resistance  data  the  hard  surface  rolling  resis¬ 
tance  and  any  zero  offset  (remaining  from  the 
static  calibration)  in  the  load  cells. 

In  tiiese  tests  the  vehicle  was  towed  with  all  four 
hubs  locked  out.  During  initial  calibration  tests. 
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Figure  1.  CRREL 
instrumented  vehicle. 


the  effects  of  the  disk  brake  pads  on  the  hard 
surface  rolling  resistance  values  were  examined. 
As  long  as  the  brakes  were  not  applied  during  the 
hard  surface  tests,  or  during  any  subsequent  tests, 
the  calibration  and  hard  surface  values  remained 
valid.  It  did  not  matter  if  the  pads  were  installed  or 
not,  although  the  hard  surface  resistance  is  higher 
with  the  pads  in  place.  By  closing  the  valves  in  the 
CIV  brake  system,  inadvertent  braking  can  be 


avoided.  Of  course,  with  no  effective  brakes  on  the 
vehicle,  care  must  be  taken  with  the  tow  vehicle 
and  tow  hitch. 

Most  tests  in  undisturbed  snow  were  conducted 
by  towing  the  vehicle  at  about  5  km/hr  using  the 
off-set  tow  hitch  shown  in  Figure  2.  A  few  tests 
were  done  without  towing  the  CIV,  using  the  front 
wheels  to  drive  the  vehicle  in  deep  snow.  Addi¬ 
tionally,  some  tests  were  conducted  with  the  CIV 


Figure  2.  Off-center 
towing  hitch. 
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Table  2.  Tire  deflection  data. 


Wheel  location 

Nominal 

load 

(N) 

Deformed 

radius 

(cm) 

Maximum 

undeformed 

radius 

(cm) 

Maximum 

section 

width 

(cm) 

Inflation  pressure  103  kPa 

left  front 

6800 

30.5 

36.5 

27.2 

right  front 

30,0 

36.5 

27.2 

left  rear 

5510 

31.5 

37.0 

27.2 

right  rear 

32.5 

37.0 

27.2 

Inflation  pressure  179  kPa 

left  front 

6800 

33.0 

36.7 

25.6 

right  front 

32.4 

37.5 

26.3 

left  rear 

5510 

33.5 

37.5 

25.5 

right  rear 

34.2 

38.0 

25.5 

following  in  the  ruts  made  by  the  tow  vehicle. 
After  each  test  run,  the  imdisturbed  snow  depths 
along  both  the  left  and  right  wheel  ruts  were 
measured  at  1-m  intervals.  These  data  were  later 
averaged  to  obtain  an  undisturbed  snow  depth 
value  for  each  side  of  the  vehicle.  At  the  conclusion 
of  a  series  of  tests,  several  snow  pits  were  used  to 
characterize  the  snow  and  to  determine  an  aver¬ 
age  density  for  each  test. 


RESULTS  AND  ANALYSIS 

Table  3  contains  a  summary  of  the  data  col¬ 
lected.  The  data  are  separated  into  sets  based  on 
the  day  of  the  test  and  the  tire  inflation  pressure. 
Initial  snow  depths  ranged  between  10  and  37  cm 
and  average  initial  snow  densities  ranged  from 
0.11  to  0.25  Mg/m^.  All  of  the  snow  data  collected 
are  given  in  Appendix  A. 


where  /i  is  the  snow  depth,  and  pf  is  the  theoretical 
final  density,  which  is  determined  from  the  fol¬ 
lowing  with  pmax  the  maximum  contact  pressure: 


pf  =  0.50  Mg/m^ 
pf  =  0.55  Mg/m^ 
pf  =  0.60  Mg/m^ 
pf  =  0.65  Mg/m^ 


for  Pmax  —  210  kPa 
for  Pmax  >  210  kPa 
for  Pmax  >  350  kPa 
for  Pmax  >  700  kPa. 


A  value  of  0.50  Mg/m^  is  used  for  the  CIV  with  the 
inflation  pressure  and  tires  used  in  these  tests. 

Three  equations  are  plotted  on  Figure  3.  These 
are 


Rs  =  15.62  (2) 


which  is  a  fit  of  the  data  collected  during  the  CIV 
tests  reported  herein 

Rs  =  11.25  (poflw)^’^  (3) 


Leading  tire  resistance 

The  leading  tire  data  from  Table  3  are  plotted  in 
Figure  3  using  the  resistance  parameter  (poaw)  in 
kilograms  per  meter,  developed  by  Ridunond  et 
al.  (1990).  The  variablesin  this  parameter  are  illus¬ 
trated  in  Figure  4;  po  is  the  average  initial  rmdis- 
turbed  density  of  the  snow,  a  is  the  arc  length  of  the 
tire  in  contact  with  the  snow,  and  w  is  the  maxi¬ 
mum  tire  width  in  meters.  The  arc  length  a  is 
determined  from  the  sinkage  z  £>nd  the  tire  radius, 
where  the  sinkage  is  determined  from 


which  was  obtained  for  tiie  CIV  alone  in  1988  and 
1989  (Richmond  et  al.  1990)  and 

Rs  =  68.083  (poflo;)”-^^^  (4) 

which  is  the  currently  used  shallow  snow  resis¬ 
tance  equation  (Richmond  et  aL  1990).  Rs  is  the 
resistance  per  wheel  (in  Newtons)  for  a  vehicle 
traveling  dvrough  primarily  undisturbed  snow. 
Combining  the  old  and  new  CIV  data  yields  the 
following  equation 


z  =  h\ 


^  n  ^ 

l_Po 


V 


pf 


(1) 


Rs  =  21.313  (poflw^)^’^^^®  (5) 

which  has  an  value  of  0.58  compared  to  0.76  for 
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Figure  3.  Leading  tire  motion  resistance. 


eq  3.  Further  analysis  of  this  data  set  (old  and  new 
CIV  data)  using  a  stepwise  regression  procedure 
showed  that  an  equation  including  only  snow 
depth  and  the  ratio  of  arc  length  «  to  1/4  the 
perimeter  of  the  wheel  provided  a  slightly  better 
correlation  to  the  data.  This  equation 


Rs  =  66.148  +  51.915/z  -  490.26 


(6) 


has  an  of  0.68.  It  is  interesting  that  neither  the 
wheel  width  nor  snow  density  is  included  in  the 
equation,  even  though  wheel  width  varied  from 
0.156  to  0.274  m  (albeit  there  are  only  a  few  data 
points  at  the  narrowest  width).  The  lack  of  a  wheel 
width  effect  was  also  suggested  by  the  analysis  of 
snow  deformation  by  a  wheel  presented  in  Ap¬ 


pendix  B.  Snow  density  (i.e.,  snow  characteristics 
other  than  depth)  may  have  been  too  similar  be¬ 
tween  tests  to  develop  a  clear  effect. 

Trailing  tire  resistance 

The  initial  analysis  consisted  of  examining  scat¬ 
ter  plots  of  the  independent  variables  (snow  depth, 
density,  tire  width,  sinkage  and  velocity)  and  sev¬ 
eral  different  combinations  of  these  parameters  vs. 
several  forms  of  the  dependent  variable:  rear 
wheel  motion  resistance  (Rj).  The  following  forms 
of  the  trailing  wheel  motion  resistance  were  con¬ 
sidered — rear  tire  motion  resistance  (Rr),  rear  tire 
motion  resistance  divided  by  vertical  load 
(coefficient  form),  ratio  of  rear  tire  motion  resis¬ 
tance  to  front  tire  motion  resistance  (R^  /Rf),  and 
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Table  4.  Average  rear  resistance  coefficient  and 
standard  deviations  by  data  set. 


Data 

set 

Number  of 
observations 

Avg.  resistance 
coefficient 

Standard 

deviation 

Average 
snow  depth 
(cm) 

Average 
snow  density 
(Mglm^) 

1 

4 

0.0207 

0.0095 

15.6 

0.165 

2 

4 

0.0171 

0.0044 

17.6 

0.160 

3 

8 

0.0056 

0.0036 

20.0 

0.220 

4 

12 

0.0140 

0.0061 

22.2 

0.220 

5 

6 

0.0728 

0.0120 

33.8 

0.160 

6 

4 

0.0068 

0.0039 

12.6 

0.110 

7 

6 

0.0199 

0.0072 

15.4 

0.110 

8 

2 

0.0152 

— 

15.8 

0.250 

9 

4 

0.0378 

0.0168 

14.6 

0.230 

the  ratio  of  rear  tire  motion  resistance  to  front  tire 
motion  resistance  in  coefficient  form  {RrVf  /RiVr)> 
This  analysis  was  done  set  by  set  and  for  all  of  the 
data  combined.  No  clear  dependence  was  found 
between  any  of  the  variables  and  any  of  the  rear 
tire  motion  resistance  "parameters." 

It  was  next  hypothesized  that  the  rear  tire  mo¬ 
tion  resistance  coefficient  (Rr  r)  ^  con¬ 

stant  (at  least  for  the  CIV)  under  similar  test  condi¬ 
tions.  To  investigate  this  hypothesis,  the  average 
and  standard  deviation  for  each  data  set  were 
calculated  (Table  4);  those  sets  that  were  collected 


under  clearly  different  conditions  were  not  in¬ 
cluded  in  the  analysis.  Data  set  5,  which  clearly  has 
a  higher  resistance  coefficient,  was  neglected  since 
it  was  obtained  in  snow  that  was  much  deeper 
than  that  in  the  other  data  sets.  Data  sets  6  and  8 
were  obtained  by  towing  the  CIV  in  the  ruts  that 
were  made  by  the  tow  vehicle,  thus  the  trailing  tire 
data  are  actually  the  fourth  pass  over  the  snow 
cover. 

Inspection  of  the  remaining  data  sets  (1, 2, 3, 4, 
7  and  9)  showed  that  the  average  values  of  the 
highest  (set  9)  and  lowest  (set  3)  didn't  seem  to 
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Table  5.  Analysis  of  variance. 

Data  Number  of 

set  observations  Mean 


1 

2 

3 

4 
7 
9 


4  0.0207 

4  0.0171 

8  0.0056 

12  0.0140 

6  0.0199 

4  0.0378 


Total  38 


For  all  data  sets 


SS 

df 

MS 

F 

TOTAL 

0.016013 

38 

CF 

0.01059 

1 

TR 

0.002986 

5 

0.000572 

RES 

0.002437 

32 

0.0000761 

7.516 

Since  F0.95 

(5,32)  =  2.5,  there  is  a  difference  between  data  sets. 

For  data  sets  1,2,3,4  and  7 

SS 

df 

MS 

F 

TOTAL 

0.009175 

34 

CF 

0.0002698 

1 

TR 

0.0075927 

4 

0.001898 

RES 

0.0013125 

29 

0.000045286 

41.92 

Since  Fo.95(4,29)  =  2.7,  there  is  a  difference  between  data  sets. 


For  data  sets  1,2,4  and  7 


SS 

df 

MS 

F 

TOTAL 

0.008817 

26 

CF 

0.007399 

1 

TR 

0.0002128 

3 

0.00007093 

RES 

0.001205 

22 

0.00005479 

1.295 

Since  Fo.95(3,22)  =  3.05,  there  is  no  difference  between  data  sets. 


agree  with  the  other  data  (Table  4).  There  was  no 
apparent  reason  for  this  on  the  basis  of  test  condi¬ 
tions  or  procedure,  although  data  set  9  may  have 
been  collected  in  hard  "crusty"  snow.  To  legiti¬ 
mately  combine  these  data  sets,  the  data  need  to  be 
random  and  normally  distributed  and  the  vari¬ 
ances  need  to  be  homogeneous.  Once  this  is  shown, 
an  analysis  of  variance  (ANOVA)  test  or  multiple 
range  test  can  be  used  to  determine  if  the  data  sets 
can  be  combined. 

Figure  5  is  a  normal  probability  plot  of  data  sets 


1, 2, 3, 4, 7  and  9.  The  data  fell  on  nearly  a  straight 
line,  indicating  that  they  are  random  and  normally 
distributed.  Bartlett's  test  for  homogeneity  of  the 
error  variance  showed  that  the  variances  were  not 
different  at  the  95%  confidence  level.  A  multiple 
range  test  indicated  that  sets  3  and  9  should  not  be 
included,  and  Table  5  contains  the  results  of  a 
series  of  ANOVA  tests  to  determine  which  of  the 
data  could  be  combined. 

Between  each  ANOVA  test,  a  procedure  fol¬ 
lowing  Grubbs  (1969)  was  used  to  determine  if  the 
extreme  values  could  be  considered  outlying  ob¬ 
servations.  This  analysis  also  showed  that  sets  3 
and  9  can  be  dropped  from  further  analysis  on  the 
basis  of  statistical  arguments.  The  average  rear  tire 
motion  resistance  coefficient  Rr/Vr  of  the  remain¬ 
ing  data  sets  is  0.0169  with  a  95%  confidence 
interval  of  0.0138  to  0,0199.  These  values  are  shown 
graphically  on  Figure  6.  It  should  be  noted  that  if 
both  sets  3  and  9  are  included,  the  average  Rr/Vr 
is  0.0167  with  a  confidence  interval  of  0.0127  to 
0.0207,  which  really  isn't  much  different  from  the 
values  obtained  from  the  four  data  sets,  although 
it  is  serendipitous  that  the  effects  of  these  two  tests 
(3  and  9)  average  each  other  out. 

From  the  above  analysis,  it  appears  that  a  trail¬ 
ing  tire  in  dry  snow  that  is  less  than  0.22  m  deep 
and  that  has  a  density  less  than  about  0.250  Mg/nr 
will  have  a  resistance  coefficient  of  about  0.017, 
although  values  as  low  as  0.0  and  as  high  as  0.64 
were  measured.  No  data  were  obtained  for  snow 

Q 

densities  greater  than  0.250  Mg/ m  .  It  should  also 
be  noted  that  the  rear  wheels  of  the  CIV  carry  less 
weight  than  the  front  wheels,  and  the  effect  of  a 
trailing  tire  carrying  a  higher  weight  is  not  known. 

Deep  snow 

Tests  0125f-0125i  (data  set  5)  were  done  in 
relatively  deep  snow  compared  to  the  other  tests. 
These  deep  snow  tests  resulted  in  limited  data  for 
two  different  trailing  tire  configurations.  Data  were 
obtained  similar  to  those  described  above  (all  four 
wheels  of  the  CIV  rolling  freely  in  undisturbed 
snow)  and  with  the  front  wheels  driving  while  the 
rear  wheels  were  rolling  free.  This  second  condi¬ 
tion,  freely  rolling  wheel  trailing  a  driven  tire,  is 
significantly  different  from  that  described  earlier. 

To  tow  the  CIV  through  this  deep  (36  cm)  snow, 
a  packed  path  was  first  made  for  the  lead  vehicle. 
Additionally,  the  snow  was  deep  enough  such 
that  the  undercarriage  of  the  lead  vehicle  dis¬ 
turbed  the  snow  in  which  the  right  side  wheels  of 
the  CIV  traveled.  Thus,  in  Table  3,  for  tests  0125f 
and  0125g,  data  are  only  presented  for  the  left  side 
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Trailing  Tire  Resistance  Coefficient,  R/V|, 


Data  Set  Number 

Figure  6.  Ninety-five  percent  confidence  intervals  for  the  means.  Shaded  area  is  the  mean 
and  95%  confidence  interval  for  combined  sets  1, 2, 4  and  7. 


Figure  7.  Snow  in  front  of  wheel. 
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Figure  8.  CIV  in  deep  snow. 


of  the  CrV.  The  snow  in  these  tests  was  about  36  cm 
deep,  which  is  about  equal  to  the  imdeflected 
radius  of  the  tires.  Figure  7  shows  how  the  snow 
piled  up  in  front  of  a  wheel,  although  it  is  not  clear 
whether  this  snow  came  from  xmder  the  vehicle  or 
was  pushed  up  from  the  area  in  front  of  the  tire 
(note  the  "bow  wave"  in  front  of  the  CIV  in  Fig.  8). 

Because  of  the  difficulty  of  towing  the  CIV  in 
snow  this  deep,  it  was  disconnected  from  the  tow 
vehicle  and  driven  in  front-wheel  drive.  It  became 
immobilized  when  the  snow  piled  up  in  front  of 
the  bumper  (Fig.  8),  so  that  tire  chains  had  to  be 
installed  on  the  front  wheels.  This,  and  slightly 
shallower  snow  in  the  succeeding  area,  allowed 
some  data  to  be  collected.  Tests  0125h  and  0125i 
yielded  trailing  tire  resistance  behind  a  driven  tire. 
The  trailing  tire  resistance  coefficients  were  0.064, 
0,057, 0.064  and  0.077  for  the  left  and  right  sides  of 
tests  0125h  and  0125i  respectively.  These  values 
and  those  of  tests  0125f  and  0125g  are  much  higher 
than  those  measured  in  the  other  tests;  this  is  most 
likely  attributable  to  snow  faUmg  back  into  the  rut 
behind  the  front  wheels  in  deep  snow,  or  in  the 
case  of  tests  0125h  and  0125i,  attributable  to  the 
disruption  of  the  packed  snow  rut  by  a  preceding 
driven  tire.  (A  preceding  driven  tire  disrupts  the 
snow  through  the  effects  of  wheel  slip;  compared 
to  the  effects  of  a  freely  rolling  wheel,  the  snow  in 
the  bottom  of  the  rut  is  made  uneven  and  looser.) 

Blaisdell  (1987)  measured  trailing  tire  resistance 


using  the  CIV  with  load  cells  mounted  only  on  the 
front  axle.  He  used  two  different  indirect  methods 
to  accomplish  this.  In  one  procedure,  he  drove  the 
front  wheels  at  a  steady  speed  and  assumed  that 
the  measured  force  on  the  front  load  cells  was 
equal  to  the  trailing  tire  resistance  (the  measured 
forces  for  the  CIV  are  analyzed  below).  A  compari¬ 
son  of  these  measurements  are  made  with  those  of 
tests  0125h  and  0125i  in  Figure  9,  which  shows  all 
the  values  falling  in  the  same  range. 

To  see  if  the  rear  coefficients  were  related  to  the 
slip  of  a  preceding  driven  wheel,  the  leading  tire's 
Differential  Interface  Velocity  (DIV)  was  compared 
with  the  trailing  tire  resistance  coefficients.  There 
did  not  appear  to  be  any  clear  cause  and  effect. 

Undercarriage  drag 

Figure  10  shows  a  force  balance  on  the  entire 
vehicle  (neglecting  air  resistance);  once  the  vehicle 
is  moving  at  a  steady  speed,  the  sum  of  forces 
should  be  zero.  If  the  CIV  is  driven  at  a  steady 
speed,  the  load  cells  on  the  driven  axles  measure 
the  resistance  to  motion  felt  by  the  rest  of  the 
vehicle  (i.e.,  the  net  traction  measured  by  these 
load  cells  must  equal  the  sum  of  all  the  other 
resistances).  Since  some  of  the  CIV  undercarriage 
was  dragging  in  the  deep  snow  tests,  undercar¬ 
riage  drag  can  be  estimated.  Minimum  ground 
clearance  for  the  CIV  is  14.0  cm  for  the  torsion  bar 
brackets  and  19.1  cm  for  the  front  differential,  and 


11 


Source:  Blaisdell  (1987)  0125h-l  0125h-r  0125i-l  0125i-r 


Figure  9.  Comparison  of  values  measured  by  Blaisdell  (1987)  and  those  obtained  here  for 
trailing  tire  resistance  behind  a  driven  wheel.  Blaisdell's  data  are  shown  ±  1  standard 
deviation. 


vehiclesmkagewasabout28cminthedeepsnow.  of  the  sum  of  the  load  cell  output,  reading  by 
Table6isasummaryoftheforcesontheCIV;hard  reading,  and  the  corresponding  velocity.  It  does 
surface  resistance  was  estimated  using  equations  not  appear  that  momentum  has  an  effect,  since  the 

developed  by  Shoop  (1992)  for  the  left  and  right  velocities  are  nearly  steady,  and  the  small  varia- 

front  tires  of  the  CIV  based  on  velocity.  Since  no  tions  in  velocity  don't  necessarily  correspond  to 

data  are  available  on  die  hard  surface  resistance  of  changes  in  the  sum  of  the  forces  on  the  vehicle. 

therearCIVtires,itwasassumedtobethesameas  The  force  ekerted  on  the  undercarriage  is  esti- 

the  front.  From  Figure  10,  the  sum  of  the  measured  mated  to  be  734  and 770 N  for  tests  0125h  and  0125i 

forces  should  equal  the  force  ascribable  to  the  respectively.  It  can  be  seen  that,  without  including 

dragging  undercarriage,  if  the  momentum  effect  theeffectof  the  slope,  the  vehicle  would  have  been 

can  be  assumed  to  be  zero  (i.e.,  the  vehicle  is  not  immobilized.  Using  the  information  from  tests 

acceleratingordecelerating).Figurellgivesplots  0215f  and  0215g  to  obtain  a  value  of  resistance 


— ►  +  Forces 
^ - Forces 


Figure  1 0.  Forces  acting  on  CIV  (front  wheel  drive).  Net  traction + ^(terrain,  internal) 
+  ^gravity)*  +  Amomentum  +  ^undercarriage)  — 

*Horizontal  component  of  gravitational  force  used  when  operating  on  a  slope. 
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Table  6.  Analysis  of  undercarriage  drag,  tests 
0125h  and  0125i. 


Test  0125h 

Test  oust 

Total  average  vertical  force 

25,371  N 

25,280  N 

Longitudinal  forces®: 

Average  left  front 

454.6 

492.7 

Average  right  front 

634.0 

771.0 

Average  left  rear 

-368.3 

-368.3 

Average  right  rear 

-347.8 

-456.4 

Hard  surface  resistance*^ 

-507.4 

-505.6 

Net  force  (available)  w/o  gravity 

-134.9^ 

-66.6 

Average  force  due  to  gravity^ 

869 

837 

Net  force  available  with  gravity® 

734 

770.4 

^From  Table  3. 

^Estimated  using  an  average  coefficient  of  0.02. 

A  negative  value  infers  that  the  vehicle  would  be  immobilized 
(no  tractive  reserve). 

^The  CIV  tilt  sensor  measured  an  average  angle  ((j))  of -1.96  and 
-1.89°,  indicating  that  the  CIV  was  moving  downhill.  The 
force  due  to  gravity  was  calculated  by  multiplying  the  total 
average  vertical  force  by  sin 

®Since  the  vehicle  was  assumed  to  be  traveling  at  a  steady 
speed,  this  value  must  be  attributed  to  undercarriage  drag. 

attributable  to  the  four  wheels  and  the  average 
undercarriage  resistance  results  in  an  undercar¬ 
riage  resistance  that  is  22%  of  the  total  vehicle 
resistance  or  1.14  times  Rs-  This  is  interesting  in 
that  Richmond  et  al.  (in  prep)  use  a  multiplier  of 
1.25  to  estimate  imdercarriage  drag  for  this  condi¬ 
tion. 

Multiple  passes 

Data  sets  6  and  8  were  obtained  by  towing  the 
CIV  such  that  the  wheels  traveled  through  the  ruts 
created  by  the  tow  vehicle.  The  forces  measured 
may  be  assumed  to  be  estimates  of  third  and 
fourth  trailing  tires.  It  was  further  assumed  that 
tests  0126g-0126i  could  represent  the  first  and 
second  wheels;  together,  these  data  represent  a 
complete  set  of  resistance  data  for  axles  1  through 
4  (Table  7).  Tests  0202f-0202h  were  analyzed  simi¬ 
larly. 

Comparing  the  average  resistance  coefficients 
between  the  two  data  sets  shows  that  the  resis¬ 
tance  data  measured  during  the  0202  tests  (snow 
depth  15.8  cm  with  a  density  of  250  kg/m^)  are 
much  higher  than  the  data  from  the  0126  tests 
(snow  depth  12.6  cm  with  a  density  of  110  kg/m^). 


This  is  most  likely  caused  by  snow  fallingback  into 
the  rut,  which  is  common  for  deep  snow  covers 
(Fig.  12)  and  may  have  happened  in  this  test. 
Unfortunately,  this  effect  was  not  identified  early 
enough  during  the  testing  to  be  fully  investigated. 
These  data  are  shown  graphically  in  Figure  13.  The 
data  of  test  0126  seem  to  agree  well  with  the  theory 
that  succeeding  tires  should  have  decreasing  resis¬ 
tance,  the  average  values  for  the  third  and  fourth 
"axles"  being  0.0094  and  0.00685  respectively. 

The  implication  of  this  analysis  is  that,  for  ve¬ 
hicles  with  more  than  two  axles  such  as  the  HEMTT 
and  LAV,  a  rationale  needs  to  be  determined  for 
modeling  these  additional  trailing  tires.  It  does  not 
appear  that  individual  wheel  resistance  on  large 
vehicles  will  be  directly  measurable  in  the  near 
future. 


SHALLOW  SNOW  RESISTANCE  MODEL 

As  mention  earlier,  Richmond  et  al.  (1990)  ana¬ 
lyzed  resistance  data  from  a  large  number  of  tests 
and  vehicles,  obtaining  eq  4  from  data  for  both 
wheeled  and  tracked  vehicles.  Richmond  (1990) 
combined  the  data  into  groups  of  like  vehicles  and 
developed  a  set  of  equations  of  the  same  form  as  eq 
4,  but  no  significant  improvement  in  the  correla¬ 
tion  coefficients  resulted.  In  these  earlier  analyses, 
only  the  data  from  the  CIV  could  be  considered  as 
leading  wheel  data.  The  data  obtained  from  the 
other  vehicles  were  "whole"  vehicle  resistances 
(from  snow)  values.  To  try  and  reduce  the  data  to 
fit  one  equation,  the  whole  vehicle  resistance  was 
simply  divided  by  the  number  of  wheels  for  each 
vehicle. 

Appendix  B  contains  a  discussion  of  the  defor¬ 
mation  of  snow  by  a  wheel  and  concludes  along 
with  eq  6  that  wheel  width  may  not  be  as  an 
important  parameter  as  arc  length  and  initial  den¬ 
sity.  Additionally,  the  ratio  of  wheel  radius  to 
sinkage  should  be  considered  to  account  for  in¬ 
creased  deformation  (more  bulldozing)  in  deeper 
snow  covers,  possibly  in  the  form  {^alm)  noted 
earlier.  With  these  ideas  in  mind,  and  with  the  new 
trailing  tire  information,  an  attempt  was  made  to 
improve  the  currently  used  algorithm  (eq  4  and  1) 
for  predicting  motion  resistance  for  wheeled  ve¬ 
hicles  in  snow. 

After  several  attempts  were  made  to  extend  the 
CrV  trailing  tire  data  to  the  HMMWV  data,  with¬ 
out  any  clear  success,  the  motion  resistance  data 
base  obtained  during  the  Wheels  vs.  Tracks  shal¬ 
low  snow  tests  was  reexamined.  Appendix  C  con- 
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Figure  11.  Net  force  and  velocity  vs.  time. 


tains  the  data  presented  by  Blaisdell  et  al.  (1990), 
Richmond  et  al.  (1990),  and  Green  and  Blaisdell 
(1991),  combined  with  the  data  presented  herein. 
The  data  are  presented  in  the  form  in  which  they 
were  used  to  develop  the  following  algorithms. 

The  first  step  was  to  reevaluate  the  sinkage 
relationship  of  Richmond  etal.  (1990).  It  was  found 
that,  if  the  tracked  vehicle  data  were  ignored,  a 
linear  trend  exists  between  the  highest  contact 
pressure  of  a  preceding  wheel  and  the  theoretical 


final  density  (pf).  The  sinkage  prediction  algo¬ 
rithm  for  undisturbed  snow  remains  as  (eq  1) 


and  the  final  density  is  now  calculated  using 

Pf  =  0.519  +  0.0023p„,,  (7) 
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Table  7.  Multiple  axle  analysis  of  resistance,  tests  0126  and  0202. 


Test 

no. 

Snow 

depth 

(cm) 

Density 

(Mg/m^) 

Farce  (N) 

Coefficient  (R/V) 

r 

2 

3 

4 

1 

2 

3 

4 

0126e-l 

12.2 

0.11 

_ 

_ 

54.3 

46.7 

_ 

— 

0.0079 

0.0088 

0126e-r 

11.2 

0.11 

— 

— 

93.0 

3.1 

- - 

— 

0.0136 

0.0006 

0126M 

14.2 

0.11 

— 

— 

74.7 

60.5 

- - 

— 

0.0107 

0.0110 

0126f-r 

12.8 

0.11 

— 

— 

104.1 

36.9 

— 

— 

0.0155 

0.0070 

0126g-l 

13.6 

0.11 

167.7 

85.0 

— 

— 

0.0241 

0.0163 

— 

— 

0126g-r 

12.7 

0.11 

182.8 

79.6 

— 

— 

0.0265 

0.0145 

— 

— 

Average 

12.7 

0.11 

175.3 

82.3 

81.5 

36.8 

0.0253 

0.0154 

0.0094 

0.0069 

0202f-l 

18.5 

0.25 

_ 

_ 

657.0 

167.3 

— 

— 

0.1002 

0.0304 

0202f-r 

13.2 

0.25 

— 

— 

311.8 

0.0 

— 

— 

0.0415 

0.0000 

0202g-l 

16.1 

0.23 

876.7 

111.2 

— 

— 

0.1271 

0.0207 

— 

— 

0202g-r 

12.8 

0.23 

521.3 

138.3 

— 

— 

0.0724 

0.0247 

— 

— 

0202h-l 

18.6 

0,23 

831.8 

327.4 

— 

— 

0.1187 

0.0638 

— 

— 

0202h-r 

10.8 

0.23 

435.0 

240.2 

— 

— 

0.604 

0.0411 

— 

— 

Average 

15 

0.23 

666.2 

204.3 

484.4 

83.7 

0.0955 

0.0376 

0.0709 

0.0152 

*Axle  number. 


where  pmax  is  the  maximum  contact  pressure  of 
the  wheels  in  kilopascals.  Figure  14  is  a  compari¬ 
son  of  the  two  prediction  methods.  Much  better 
agreement  is  obtained  using  eq  1  and  7  as  com¬ 
pared  to  eq  1  and  the  table  of  values  for  pf. 


Using  the  resistance  data  in  Appendix  C  and  a 
nonlinear  regression  analysis  combined  with  the 
assumptions  that  1)  the  resistance  of  the  first  wheel 
through  undisturbed  snow  is  a  function  similar  to 
eq4  and  2)  that  the  resistance  of  trailing  tires  can  be 


Figure  12.  Loose  snow  in  rut  in  front  of  trailing  tire. 
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represented  as  a  coefficient  multiplied  by  the  ve¬ 
hicle  wheel  load,  and  3)  using  the  sinkage  predic¬ 
tion  algorithm  of  eq  1  and  7,  resulted  in  the  follow¬ 
ing  equation  for  1/2  the  total  vehicle  resistance 
(each  axle  is  considered  to  be  symmetrical). 

Rs  =  13.604l(poa;fl)^'^^  + 

(8) 

0.0475 /irL2  +  0.0246  lirLs 


where 

ftr  =  1.0if/i<  0.22  m 

hr  =  —iih>0.22m 

"  0.22 

where 

h  =  snow  depth  (m) 

L2  =  1/2  the  load  (N)  on  the  second  axle 
L3  =  1/2  the  load  on  the  third  axle. 


Figure  14.  Comparison  of  two  sinkage  prediction  algorithms. 
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Figure  15,  Comparison  of  predicted  and  measured  resistance  values. 


The  term  (po^)  is  in  kilograms  per  meter.  It  is 
interesting  to  note  that  adding  a  term  for  the  fourth 
axle  adds  little  to  the  accuracy  of  the  equation  and 
forces  the  coefficient  for  the  third  axle  to  be  nega¬ 
tive,  a  physical  impossibiUty .  The  rationale  for  the 
ratio  of  snow  depth  to  the  depth  of  0.22  m  is  to 
account  for  the  effect  of  snow  faUingback  into  a  rut 
behind  a  leading  wheel  in  deep  snow.  The  trailing 
tire  coefficients  obtained  are  higher  than  antici¬ 
pated,  and  this  is  thought  to  be  attributable  to 
higher  loads  on  the  trailing  axles  for  most  of  the 
vehicles  in  the  data  base,  as  compared  to  the  CIV, 
which  has  less  weight  on  the  trailing  axle. 

Other  forms  of  the  term  {poV)a)  were  investi¬ 
gated  using  the  ideas  presented  earlier;  little  or  no 
improvement  w’^as  fotmd  in  the  regression  analy¬ 
sis.  Figure  15  shows  a  comparison  of  the  predicted 
vs.  measured  resistance.  It  can  clearly  be  seen  that 
eq  8  produces  excellent  results  and  a  marked  im¬ 
provement  over  earher  algorithms.  Table  Cl  con¬ 
tains  the  actual  predictions  and  residuals. 

SUMMARY 

Several  aspects  of  the  interaction  between  snow 
and  a  wheeled  vehicle  were  investigated.  The 


resistance  of  a  trailing  wheel  can  be  considered  a 
constant  in  most  cases;  however,  snow  faUingback 
into  the  rut  behind  a  leading  wheel  results  in 
higher  values  and  seems  to  be  related  to  initial 
snow  depth.  The  true  effect  of  wheel  width  on 
resistance  of  a  leading  wheel  is  stiU  not  known  and 
should  be  studied  using  a  numerical  technique 
where  other  parameters  can  be  held  constant.  Lim¬ 
ited  data  were  presented  for  undercarriage  drag, 
and  since  this  will  be  the  primary  cause  of  immo¬ 
bilization  on  level  surfaces,  it  deserves  further 
investigation. 

A  new  and  significantly  improved  algorithm 
for  predicting  sinkage  and  motion  resistance  in 
imdisturbed  snow  was  presented.  This  improved 
algorithm  should  be  implemented  in  current  cold 
regions  mobility  models. 
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APPENDIX  A:  SNOW  DATA 


Tests  0119H,  0119  0 

Density  (Mg/m Air  Temp.:  “5°C 


Ave.  Density:  0.18 


Tests  0119  G,  0119  P 


Ave.  Density:  0.16 

a.  19  January  1993. 
Figure  Al.  Snow  data. 
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APPENDIX  B:  OBSERVATIONS  OF  SNOW  DEFORMATION  BY  A  WHEEL 

The  purpose  of  this  study  was  to  examine  the  deformation  of  snow  around  a 
rolling  wheel  and  to  identify  mechanisms  and  parameters  that  could  lead  to  an 
improved  motion  resistance  model  for  shallow  snow. 

Procedure 

Deformation  experiments  were  conducted  in  areas  of  undisturbed  snow  that 
were  marked  by  chalk  dust.  Holes  in  the  snow,  which  were  punched  both  perpen¬ 
dicular  to  and  in  the  direction  of  vehicle  travel,  were  coated  with  carpenter's  chalk 
dust.  This  was  done  by  giving  a  plastic  bottle  of  chalk  dust  with  a  tapered  nozzle  a 
quick  squeeze  into  the  hole.  A  12.5-mm  (1  /2-tn.)  diameter  dowel  was  used  to  pimch 
the  holes,  which  were  placed  76,2  mm  (3  in.)  apart  (Fig.  Bl).  After  emplacing  the 
dust,  a  vehicle  was  driven  into  the  marked  area  such  that  the  axle  passed  through 
the  row  of  perpendicular  holes,  then  the  vehicle  was  backed  out.  Whenever  possible, 
the  wheel  was  undriven.  The  snow  was  then  carefully  removed  to  reveal  a  cross 
section  of  the  deformed  area  and  the  lines  formed  by  the  chalk  (Fig.  B2).  Density 
measurements  were  made  using  a  lOO-cm^  sampler,  which  had  a  6-  by  3-cm  cross 
section.  Air  and  snow  temperatures  were  recorded  and,  in  some  cases,  so  was  the 
grain  size. 

Three  different  vehicles  were  used,  although  most  of  the  experiments  used  the 
CRREL  instrumented  vehicle.  Tire  data  for  each  vehicle  are  in  Table  Bl. 

Results  and  observations 

Photographs  and  documentary  data  are  presented  in  Figures  B3-B5.  From  these 
figures  it  can  be  seen  that,  as  the  wheel  rolls  forward  in  the  snow  pack,  deformation 
occurs  in  all  three  directions.  In  the  following  discussion,  forward  is  the  direction  in 
which  the  vehicle  is  moving,  perpendicular  refers  to  the  direction  perpendicular  to 
the  direction  that  the  vehicle  is  moving  but  parallel  to  the  ground,  and  downward 


Figure  Bl.  Snow  marked  with  chalk-dust-filled  holes. 
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Figure  B2.  Removing  snow  to  show  deformed  chalk  holes. 


or  lower  infers  motion  towards  the  groimd.  The  deformation  around  a  rolling  wheel 
is  extremely  complex  and  there  seem  to  be  four  areas  to  consider.  These  are  the  area 
imder  the  wheel  (the  rut  area),  the  area  in  front  of  the  wheel,  and  the  upper  and  lower 
areas  adjacent  to  the  wheel  sides.  The  deformation  in  front  of  the  wheel  and  to  the 
sides  of  the  wheel  in  the  upper  part  of  the  snow  cover  are  related  to  each  other,  to 
the  snow  depth  and  the  wheel  sinkage. 

Looking  first  at  the  snow  imder  the  wheel  rut,  one  can  see  that  the  snow  at  the 
centerline  of  the  wheel  is  pushed  forward  and  downward  as  the  wheel  rolls  forward. 
Snow  that  is  off  the  centerline  is  pushed  off  at  an  increasing  angle  by  the  compacting 
snow  closer  to  the  centerline.  This  is  seen  in  the  chalk  exposed  in  the  area  of  the  wheel 
rut.  The  chevron  shapes  are  created  because  the  snow  near  the  center  of  the  snow 
pack  is  pushed  out  further  (perpendicularly  and  forward)  than  the  snow  near  the 
snow  or  ground  surfaces.  As  forward  motion  continues,  the  upper  snow  is  pushed 
downward  (see  Fig.  B3a,  c  and  B4). 


Table  Bl.  Vehicle  and  wheel  characteristics. 

CKREL  Instrumented  Vehicle  (CtV),  1977  Jeep  Cherokee 
Tires:  Michelin  LT  235XCH4, 179  kPa 
Radius:  37.5  cm 

Maximum  deformed  width:  26  cm 

CRREL  Support  Truck,  1989  Chevrolet  3/4  ton  pickup 
Tires:  LT  245/75R16  M&S 
Radius:  37  cm 

Maximum  deformed  width:  not  recorded 

HEMTT,  M997  Heavy  Expanded  Mobility  Tactical  Truck 
Tires:  Michelin  16.0  R20 
Radius:  not  recorded 

Maximum  deformed  width:  48  cm  _ 
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*Snow  density  (Mg/m^)  typical. 


5  cm 


b.  27  January  1992.  Snow  depth:  21.6  cm;  snow  temperature: -3°C; grain  size:  l.Omnfi; 
air  temperature:  -7°C;  wheel  size:  LT  235XCH4, 179  kPa;  wheel  sinkage:  16  cm. 

Figure  B3  (cont'd).  Snow  deformation  and  density  distributions. 
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c.  29  January  1992,  Snow  depth:  11,5  cm;  snow  temperature:  -0°C;  air  temperature:  -2°C;  wheel  size:  LT 235XCH4, 
179  kPa;  wheel  sinkage:  8  cm. 


d.  3  February  1992,  Snow  depth:  31,5  cm;  snow  temperature:  -l^C;  air  temperature:  1°C;  wheel  size:  LT 235XCH4, 
179  kPa;  snow  was  removed  from  the  surface  so  that  the  wheel  could  be  driven  into  the  deep  snow.  The  axle  did  not 
pass  the  chalk  dust  holes  (4-5  cm  short). 

Figure  B3  (cont'd). 
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/.  7  February  1992.  Snow  depth:  16  cm;  snow  temperature:  -1°C;  air  temperature:  -11°C;  wheel  size:  LT  235XCH4, 
179  kPa;  wheel  sinkage:  10  cm. 

Figure  B3  (cont'd). 
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g.  8  February  1992.  Snow  depth:  25  cm;  snow  temperature:  -13°C;  air  temperature:  -irC;  wheel  size:  LT  245/ 
75  R16,  M&S;  wheel  sinkage:  about  14  cm. 


Figure  B3  (coni' d).  Snow  deformation  and  density  distributions. 
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Figure  B4.  Deformation  in  direction  of  vehicle  travel,  same  snow  properties  as  in  Figure 
B3f 


For  the  snow  located  to  the  side  of  the  wheel,  it  appears  that  the  forward 
displacement  component  is  smaller  than  the  perpendicular  displacement,  since  the 
chalk  at  further  distances  from  the  centerline  has  little  forward  displacement 
compared  to  chalk  near  the  wheel  (Fig.  B3d).  This  may  occur  as  the  full  load  on  the 
wheel  is  transferred  to  the  snow,  causing  additional  plastic  deformation  in  the 
direction  perpendicular  to  the  wheel.  Maximum  deformation  seems  to  occur  at  or 
slightly  below  the  bottom  of  the  rut.  Additionally,  there  appears  to  be  little 
displacement  of  the  snow  near  the  surface  compared  to  snow  deeper  in  the  snow 
cover.  This  suggests  that  the  strength  of  the  snow  lower  in  the  snowpack  is  more 
important  in  resistance  calculations. 

As  wheel  sinkage  approaches  the  radius  of  the  wheel,  the  forward  part  of  the 
wheel  acts  more  as  a  bulldozer  than  a  compactor  (as  in  the  case  of  relatively  shallow 
snow)  and  deformation  (both  forward  and  perpendicular)  is  higher  in  the  snow  pack 
(compare  Fig.  B3b  and  c).  This  is  also  observed  in  Figures  B4  and  B5  for  snow 
deformation  in  front  of  the  wheel.  In  shallow  snow  the  deformation  area  does  not 
extend  beyond  the  leading  edge  of  the  tire  (Fig.  B5),  while  for  deeper  snow  it  does 
(Fig.  B4).  This  effect  can  be  illustrated  by  examining  the  velocity  vector  components 
of  a  rigid,  non-slipping  wheel.  Figure  B6  shows  the  forward,  vertical  and  resultant 
vectors.  One  can  see  how  the  forward  component  is  reduced  as  one  moves  along  the 
tire  surface  towards  the  supporting  surface. 

Most  of  the  experiments  were  done  in  relatively  low-density  snow.  When  the 
snow  had  a  higher  density,  greater  deformation  was  observed  (compare  Fig.  B3b 
and  c).  This  agrees  with  observations  of  soil  compaction  and  suggests  that  there  may 
be  a  controlling  density  (maximum  packing  density  without  significant  crystal 
deformation).  Recent  analysis  by  Shapiro  et  al.  (in  prep.)  suggests  that  intergranular 
bonding  (between  snow  crystals)  is  more  important  than  previously  suspected,  and 
will  greatly  control  this  type  of  displacement.  There  is  currently  no  proven  method 
of  measuring  bond  strength  and  density  remains  the  primary  field  index  of  snow 
strength. 
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Fi<mre  B5  Deformation  (horizontal  and  lateral)  and  initial  snow  density,  8  February  1992.  Snow  depth:  14.5  cm; 
snoii;  temperature:  -13°C;  air  temperature:  -11°C;  wheel  size:  LT  245/75  R16,  M&S;  wheel  sinkage:  10.5  cm. 
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a.  Forward  component  unit  vectors. 


b.  Vertical  component  unit  vectors. 


Figure  B6.  Forces  on  a  rigid  rolling  wheel 


Two  tests  were  done  in  fairly  deep  snow  with  significantly  different  wheels. 
Figures  B3c  and  d  show  the  results  for  the  CIV  wheel  (tire  width  of  26  cm)  and  the 
HEMTT  (tire  width  of  48  cm).  For  both  tests  the  maximxim  perpendicular  displace¬ 
ment  is  about  18-20  cm.  This  indicates  that  the  effect  of  tire  width  on  resistance  may 
not  be  linear. 

Summary 

The  results  of  some  simple  tests  to  examine  the  deformation  of  snow  aroimd  a 
rolling  wheel  were  presented.  Although  they  were  primarily  qualitative,  some 
interesting  phenomena  were  observed.  Snow  deformation  caused  by  a  tire  is  three- 
dimensional  and,  depending  on  the  depth  and  density,  will  primarily  occur  at  or 
below  the  wheel  rut  in  shallow  snow  and  in  the  upper  area  if  the  snow  is  deep. 
Increasing  wheel  width  does  not  seem  to  increase  deformation  in  the  direction 
perpendicular  to  vehicle  motion,  at  least  in  the  one  comparison  available  from  these 
tests.  The  ratio  of  wheel  radius  to  sinkage  should  be  considered  to  account  for 
increased  deformation  in  deeper  snow  covers. 
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APPENDIX  C:  WHEELED  VEHICLE  MOTION  RESISTANCE  DATA 


Table  Cl  contains  the  data  for  the  CIV  reported  here  plus  the  wheeled  vehicle  data 
presented  in  BlaisdeU  et  al.  (1990),  Richmond  et  al.  (1990)  and  Green  and  BlaisdeU 
(1991).  Some  corrections  were  made  to  the  data  presented  in  the  earlier  reports.  The 
resistance  values  are  presented  as  the  measured  value  minus  the  hard  surface 
resistance  divided  by  two,  and  thus  represent  the  1/2  of  the  resistance  attributable 
only  to  the  snow.  For  each  vehicle,  normal  loads  are  1/2  the  axle  loads;  some  of  the 
CIV  resistance  data  are  for  the  first  axle  only  and  thus  it  is  considered  a  one-axle 
vehicle.  The  sinkage  is  calculated  using  eq  1  and  7.  The  shaded  data  were  not 
included  in  the  regression  analysis  for  the  reasons  stated  in  the  Table  Cl  notes. 
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Table  C1.  Wheeled  vehicle  motion  resistance  data. 
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